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Abstract 



^ , We derive equation describing distribution of energy losses of the 

Q I particle propagating in fractal medium with quenched and dynamic 

, ^, ' heterogeneities. We show that in the case of the medium with fractal 

dimension 2 < D < 3 the losses A are characterized by the sublinear 
_ anomalous dependence A ~ x" with power-law dependence on the 

^fs ■ distance x from the surface and exponent a = D — 2. 

oo ■ 

^ I This letter is devoted to studying statistical properties of the coUisional 

energy losses suffered by a high energy particle passing through a randomly 

Q . inhomogeneous disordered medium. Exploration of the properties of complex 

media with the help of test particles propagating through it is one of the 
most important scientific instruments used in physics allowing, in particular, 
to study a response of the medium to particle beams or radiation coming 

^ ■ through the medium under study. A range of possible applications is quite 

j^ ■ broad, from the physics of high energy collisions to polymer physics, to name 

the few. In all cases the main quantity under study is the distribution of 
energy losses of the test particle studied as a function of the distance covered 
by the test particle in the medium. 

The original context of the problem was related to examination of ion- 
ization losses of high energy particles in ordinary homogeneous amorphous 
matter (energy straggling), see e.g. [I]. The microscopic picture underlying 
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the energy losses in this case was, evidently, that of a series of inelastic colli- 
sions of the projectile with atomic electrons resulting in excitation/ionization 
of the corresponding atoms. In between the scattering events the projectile 
trajectory is ballistic and its energy does not change. In the continuum limit 
the problem can be reformulated, for a high energy projectile and energy 
losses small compared to the energy of the incident particle, in terms of a 
one-dimensional kinetic equation suggested by Landau [2], in which the role 
of time is played by the distance along the straight line trajectory covered 
by the projectile in the medium. 

The focus of the present paper is on the statistical properties of energy 
straggling in matter characterized by strong random inhomogeneities with 
the special emphasis on the case of fractal medium. Our consideration of the 
microscopic energy loss model in a strongly inhomogeneous random medium 
with power-like correlations is realized at the level analogous to the Contin- 
uous Random Walk (CTRW) [3] that was, in particular, actively explored to 
describe thermal diffusion of particles in a fractal medium. |1H7] Physically, 
the main difference between the stochastic process describing the energy loss 
and the CTRW is in straight line trajectory of high energy particle, while 
CTRW describes haotic trajectories of particle induced by random Brownian 
forces. From mathematical point of view the distinction of these processes 
is that the energy loss problem corresponds to a sum of random positive 
quantities - energy losses at scattering events, whereas random particle dis- 
placements in CTRW can have any sign. 

Let us formulate the basic microscopic model of energy loss studied in this 
paper and consider a high energy particle incident on the medium containing 
randomly placed scattering centers. In the high energy approximation the 
trajectory of the projectile is a straight line and the particle is assumed to 
interact with all the scattering centers that happen to lie at the projectile's 
trajectory. More precisely, let us consider the particle with high energy Eq 
entering the medium at the point x = 0. Our main goal is to compute the 
distribution / (A, x) of its energy loss A at some depth x from the surface. 

The energy loss A at the point x in a given event is fully characterized by 
the set of energy losses at each of scattering event {Aj}, i = 1, ■ ■ ■ , n so that 
A = Ai -|- A2 + ■ ■ ■ -|- A„. The corresponding configuration of the scattering 
centers is, in turn, fully specified by the set of distances h,h,--- JnJn+i, 
where h is the distance between the surface and the first scattering cen- 
ter, {li} are the distances between the i'th and i + I'th scattering centers 
and, finally, In+i is the distance between the last scattering center and the 



observation point x {x = li + I2 + ■ ■ ■ + ^n+i), see Fig. [H 
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Figure 1: Particle propagating in the medium looses the energy A = Ai + 
A2 + ■ ■ ■ + A„ in n scattering events taking place at points separated by 
intervals {/j} along the trajectory of the projectile. 

In the assumption of independent energy losses at different scattering 
centers with the distribution of energy losses w{e) the probability distribution 
/(A,x) of the cumulative energy loss A at some point x is fully described 
by the probability densities ip{l) of having a spatial distance I between two 
scattering centers, see Fig. [H For large number of scattering events the energy 
losses do not depend on the probability density ^{li) of the first event, and 
in follows we will simply assume that ^{li) = 'ipih)- The problem of random 
loss is similar to that of the usual random walk, with the distance x and local 
energy loss A replacing the time and displacement correspondingly. From 
this analogy the energy loss distribution function / (A,x) is described by 
CTRW-like equation: 

px /-A 

/ (A,x) = (5 (A) * (x) + / dx'tp {x - x') / dew (e) / (A - £,x') , (1) 
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where \l/ (x) = J dytply) is the probability of having no scattering events in 
the interval [0,x]. Eq. ([T]) is conveniently solved by using the double Laplace 
transform 
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fip,q)^ / rfAe-P^ / dxe-^^f{A,x) 
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Introducing the finction g (q) = ip (g) /[I — ip (g)] the corresponding equation 
for / {p, q) can be written in the form 

/ (p, q) = i/q+~9 (q) [w (p) - i] / {p, g) (2) 

which, in turn, corresponds to the following version of the original kinetic 
equation ([1]): 

PX POO 

f (A,a;) = 5 (A) + / dx'g {x - x') / dew {e) [/ (A - e, x') - f (A, x')] (3) 
Jo Jo 



The function g (r) is found by inverse Laplace transform of the function 
g {q) and it has the meaning of the average density of scatterings along the 
direction of particle propagation at the distance r from the last scattering. 

In general, the function g (r) depends on characteristics of the medium, 
and it can be related to the so-called structure function of the medium 

G ir) = {J2^^,H^^ - ^^+n - r)) , (4) 

where Xj are coordinates of the i-th scattering center, by 

g (r) = a'G (r) , (5) 

a^ is the scattering area of the particle. In many important cases of the 
scattering of particles in a complex heterogeneous medium a microstructure 
of the medium remains unknown, but its structure function G (r) can be 
directly measured experimentally. In these cases the kinetic equation ([3]) 
can be used to predict the spectrum of energy losses in such a medium. 
The knowledge of the spectrum is extremely important, for example, in the 
problem of radiation damage of the medium which is determined not only 
by total adsorbed energy but also by the shape of the distribution of energy 
losses. 

Eq. ([3]) may be considered as a generalization of the Landau equation for 
ionization losses in amorphous media [1], which can be written in integral 
form: 

f (A,x) = 6iA) + - / dx' dew (e) [f (A - e, x') - f (A, x')] (6) 
a Jo Jo 

This equation is obtained from Eq. ([3]) in the case of the constant linear 
density of scattering centers g {x — x') = 1/a. 

Analytical solution of generalized kinetic equation ([3]) can be found in 
the case of scattering medium with fractal dimension D, when the Laplace 
transform of the function g (r) has the form: 

~giq) = iaq)-\ a = D-2 (7) 

In the infrared limit x ^ a one can approximate the Laplacw transform of 
the function w (e) hj w (p) ~ 1 — pe, where e = J ew{e)de is the average 
energy loss in a scattering event. Calculating the inverse Laplace transform 
of Eq. ([2]), we get 
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where W^ is the Wright function [S] 
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Using the distribution function (|H]) one can compute the average energy loss 
at some depth x: 

(A (x)) = [e/T (1 + a)] {x/aY , < a < 1 (10) 

We conclude that in the case of the fractal medium the average energy loss 
is characterized by fractional sublinear dependence on the distance. The 
equations (IHfTO!) constitute the main result of the paper. 

Studied in this paper problem of random energy loss can be considered 
as "dual" to the problem of thermal random walk motion of low energy 
particles. The thermal diffusion of particles in a fractal medium is known to 
be anomalously slow, see e.g. [9]. The reason of the subdifFusional motion in 
fractal porous medium is that a particle is trapped in dead end pores and 
bottlenecks, so that diffusion is slowed down and becomes anomalous. The 
physics of anomalous random energy loss in the fractal medium is different 
and is related with the presence of long-range correlations in positions of 
scattering centers. Approach developed in this work can also be apphed 
to describe particle propagation in the system with dynamic heterogeneities 
formed at the critical point of phase transition. 
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